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Abstract 

We study the evolution of homogeneous and isotropic, fiat cosmological models within 
the general scalar-tensor theory of gravity with arbitrary coupling function and potential. 
After introducing the limit of general relativity we describe the details of the phase space 
geometry. Using the methods of dynamical systems for the decoupled equation of the 
Jordan frame scalar field we find the fixed points of flows in two cases: potential domina- 
tion and matter domination. We present the conditions on the mathematical form of the 
coupling function and potential which determine the nature of the fixed points (attractor 
or other). There are two types of fixed points, both are characterized by cosmological 
evolution mimicking general relativity, but only one of the types is compatible with the 
Solar System PPN constraints. The phase space structure should also carry over to the 
Einstein frame as long as the transformation between the frames is regular which however 
is not the case for the latter (PPN compatible) fixed point. 
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1 Introduction 



Scalar-tensor theories (STT) of gravitation [TJ [2l E] have emerged in different contexts of 
theoretical physics, e.g. in Kaluza-Klein type unifications, supergravity, and low energy ap- 
proximations of string theories. The scalar-tensor action functional has been also used in 
topology as a mathematical tool for solving the geometrization conjecture of 3-manifolds via 
the Ricci flow [J]. 

In cosmology STT has been invoked to model the accelerated expansion of inflation [5] and 
dark energy [6] (for some recent papers see e.g. Ref. [7j). However, observations in the Solar 
System tend to indicate that in the intermediate-range distances the present Universe around 
us is successfully described by the Einstein tensorial gravity alone [8]. This means that only 
such models of scalar-tensor gravity are viable which in their late time cosmological evolution 
imply local observational consequences very close to those of Einstein's general relativity (GR) 

The importance of the convergence of the STT solutions to those of GR at late times, guar- 
anteeing agreement with the present post-Newtonian observations was recognized by Beken- 
stein and Meisels \±0\ already long ago. Damour and Nordtvedt [TU [12] argued that for a wide 
class of homogeneous and isotropic scalar-tensor cosmological models there exists an attrac- 
tor mechanism taking them to the limit of the corresponding general relativistic models. In 
subsequent papers the issue has been pursued by various approaches [T3J [HI [151 EH] • 

The methods of dynamical systems provide natural tools to analyze the problem. Most 
previous studies which have considered STT cosmology as a dynamical system have focused 
upon examples with specific coupling functions, see Refs. [171 HH [191 120] for works based on 
the Einstein frame and Refs. [2T1 [221 l23"l I2H [20] for the Jordan frame. In the Jordan frame the 
main properties of the general phase space geometry were outlined by Faraoni [25] , while Refs. 
|26[ [27] consider general Jordan frame dynamics with particular attention to the de Sitter fixed 
point. 

Although the equations of motion in the Einstein frame may have a simpler form, the cos- 
mological observations are easier to interpret in the Jordan frame where matter is minimally 
coupled. In general there is a one-to-one correspondence between the phase spaces in the 
two conformal frames [281 EHJ El [301 EI], provided that the transformations from one frame 
to another are regular. However, as it was pointed out in Ref. [20], the differential equation 
relating the scalar fields in these two frames becomes singular at the limit of general relativ- 
ity. Therefore utmost care must be exercised when addressing the question whether general 
relativity is an attractor, as the answer may depend on the frame chosen. 

In this paper we take the Jordan frame and consider general scalar-tensor theories which 
contain two functional degrees of freedom, the coupling function and the scalar potential 
V(^/). We perform the dynamical systems analysis for the flat Friedmann-Lemaitre-Robertson- 
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Walker (FLRW) backgrounds with ideal barotropic fluid matter. Our strategy is to find the 
fixed points for the scalar field dynamics and compare these with the conditions of the limit 
of general relativity in the Solar System, as established by the parameterized post-Newtonian 
(PPN) formalism. Therefore, if the functional forms of u(^) and V(ty) are specified from 
some considerations (e.g. the compactification manifold), our results allow to determine the 
fixed points along with their type and thus immediately decide whether general relativity is 
an attractor, i.e. whether the model at hand is viable or not. 

The plan of the paper is the following. In section 2 we write the field equations as a 
dynamical system and find the conditions which reduce these equations to the FLRW equations 
of general relativity (possibly with a cosmological constant), concluding that these conditions 
are marginally stricter than the standard PPN conditions. In section 3 we describe the general 
phase space and its physical subspace determined by the Friedmann constraint equation. In 
section 4 we consider the case of potential domination with (effectively) vanishing matter 
density, find the fixed points and determine their properties. For investigating the case of 
matter domination (with vanishing scalar potential) in section 5, we follow Refs. [12] and [15] 
to introduce a new time parameter which allows us derive a decoupled "master" equation for 
the scalar field and again perform the fixed point analysis. In section 5 we briefly discuss 
comparison to some previous works and end with a speculation about the possible relevance of 
our results for the moduli stabilization problem in string theory. Section 6 gives a summary. 

2 Scalar-tensor cosmology as a dynamical system and 
the limit of general relativity 

We consider a general scalar-tensor theory in the Jordan frame given by the action functional 

)• (1) 

Here u(^>) is a coupling function and V(ty) is a scalar potential, denotes the covariant 
derivative with respect to the metric g^, k 2 = 8tiGn, and S m is the matter part of the action 
as all other fields are included in \ m . In order to keep the effective Newtonian gravitational 
constant positive [3U [H] we assume that < \& < oo. 

The field equations for the flat (k = 0) FLRW line element and perfect barotropic fluid 
matter, p = wp, read 

H = -^ + 6^ w + yf + y4^' (2) 



S = 2^ 1 d ' x ^ 



VR{g) - ^^V P ^V P ^ - 2k 2 V(*) 
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H = a/a. The matter conservation law is the usual 

p + 3if (w+ 1) p = 



2w(tf) + 3 



'1 — 3w) p 



(4) 



(5) 



and we assume p > 0. Eqs. ((21)-© are too cumbersome to be solved analytically, but useful 
information about the general characteristics of solutions can be obtained by rewriting ©-© 
in the form of a dynamical system and finding the fixed points which describe the asymptotic 
behaviour of solutions. 

The phase space of the system is spanned by four variables H, p}. Defining ^ = 

x, ^ = y the dynamical system corresponding to equations ©-© can be written as follows: 
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(8) 
(9) 



Based on these equations we may make a couple of quick qualitative observations about 
some general features of the solutions. For example, the limit ^ — > in general implies 
\H\ — ► oo, hence the solutions can not safely pass from positive to negative values of \& (from 
"attractive" to "repulsive" gravity), but hit a space-time singularity as the curvature invariants 
diverge. Similarly, the limit 2a; +3 — ► implies \H\ — > oo with the same conclusion that passing 
through o;( 1 i r ) = — | (corresponding to the change of the sign of the scalar field kinetic term 
in the Einstein frame action) would entail a space-time singularity and is impossible. (Let us 
remark here, that these observations are quite general and do not preclude specially fine-tuned 
solutions in some fine-tuned models which may remain regular while crossing these points, like 
in Ref. E2J.) 



The limit 
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deserves a more closer examination. Let us define x+ by 

1 



2a; (x*) + 3 

Expressing H from the Friedmann constraint ©, 
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makes clear that \H\ diverg unless also y — > at the same time. What happens in 

the latter case is determined by the first term under the square root above. We can compute 
its limit by Taylor expanding 
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where Ax = rcosO, Ay = rsin6* was taken. (We have neglected the unphysical direction 



\9\ = | that corresponds to approaching the point 
where \H\ is divergent.) So, in the limit x 



0) along the line x 



x*,y 



lowest non-zero derivative of 



the value of H is determined by the 
, . If the both the first and second derivative vanish, then 

2uj(x)+3 ' 

| if | diverges implying a spacetime singularity. If the first derivative vanishes but the second 
derivative is not zero, 



2u(x) + 3 







dx 2 \2uj(x) + 3 
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then H is finite but (possibly) different for each solution as it depends on the angle of approach 
9, while the Friedmann equation in this case acquires an extra term when compared to general 
relativity. If the first derivative is not zero, 
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2uj(x) +3 
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then H approaches the value if 2 = &-(p + V(x ic )), mimicking the Friedmann equation of 
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general relativity with 8nG = — and A = — V(x+). 

To summarize, we have just observed that in the limit (a) 2 u(x)+3 ~ > ^> O 3 ) V ~ * ^ ^ ne 

^0. 



Mi)T3 °> ( b ) V 

Friedmann constraint ([TT]) tends to the form of general relativity if (c^ ,u 



(2uj(x ir )+3) 2 dx 

It must be also emphasized here that the process of taking the Taylor expansion ffl2|) hinges 
on the assumption that (d) 2 u(l)+3 * s differentiable (derivatives do not diverge) at x+. In this 
context one may also ask when the full set of equations (IHJ)- (|9j) attains the form of general 
relativity. It is easy to see that besides (a)-(d) one must also impose 
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but the latter is automatically satisfied if (c) holds, due to (1121) . (j!4p . Therefore we may 
tentatively call the conditions (a)-(d) 'the general relativity limit of scalar-tensor flat FLRW 
cosmology'. 

It is interesting to compare the cosmological GR limit to the GR limit obtained from PPN, 
which characterizes the slow motion approximation in a centrally symmetric gravitational field. 
Although the mathematical assumptions underlying the PPN formalism are clearly different 
from our cosmological reasoning above, we may still ask whether the results of both schemes 
agree with each other. In the context of PPN it is well established that the solutions of 
scalar-tensor theory approach those of general relativity when 



1 1 du 

2u{x) + 3 {2u{x) + 3) 3 dx~ (16) 

Comparison shows that the cosmological conditions (a)-(d) are marginally stricter than the 
PPN condition ( Tl6l) . since (a), (c), (d) imply that (fTBl) is satisfied, but ( TT6i) does not necessarily 
guarantee that (c) or (d) holds. 

Let us also note that there is also another special realized at 

dV(x) 

p = 0, y = 0, 2V(x.)-x. =0, (17) 

ax x » 

when the cosmological equations ©-© relax to those of general relativity featuring de Sitter 
evolution. However, as the value of ui(x t ) is not fixed by the condition (JTTJ) , this case does 
not conform with the GR limit of PPN. Therefore, even when the limits (a)-(d) and ffTTl) can 
be cosmologically indistinguishable, Solar System observations in the PPN framework can in 
principle reveal which of the two is actually realized. (In this paper when using the phrase 
'the GR limit of STT' we mean the conditions (a)-(d), as these take the STT cosmological 
equations to those of general relativity and also guarantee that the PPN condition is satisfied. 
But note that some authors, e.g. Refs. [18j[T9] have not necessarily used the same definition.) 

The general relativity limit of STT is purely given in terms of x and y. In the following 
we extract from the full dynamical system (EJ)— (EI) an independent subsystem for {x,y}, find 
its fixed points and check whether the limit of general relativity matches to an attractive fixed 
point. But before, some general remarks about the full phase space are in order. 

3 Phase space 

In the four phase space dimensions = x, 4f = y,H,p} the physical trajectories (orbits of 
solutions) are those which satisfy the Friedmann constraint ([2]), i.e. which lie on the 3-surface 



F:F(x,y,H,p)=H + H-- — u(x)- — 3^ = ■ 
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It can be readily checked that the trajectories' tangent vector, T % = (x,y, H, p), given by 
(E])-®, is perpendicular to the normal of the Friedmann surface, 



V 4 F-T^ = 0, (19) 

and therefore the system automatically obeys the Friedmann constraint, as the trajectories on 
the surface T never leave it. 

In principle the geometry of the Friedmann 3-surface in the 4-dimensional phase space is 
rather complicated to visualize, but a few general characteristics can still be given. We may 
write Eq. (fl8|) as 
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y \ 

2 X ) y ^ (2Q) 



K 2 ( P +V) iK 2 x{p+V) ' 

3x 2lu+3 

which for fixed p and x can be recognized as describing familiar conic sections: 1) for p+V > 0, 
2u + 3 > a hyperbola on the (H + j^,y) plane, 2) for p + V > 0, 2u + 3 < an ellipse also 
on the (H + y) plane, while 3) for p + V < 0, 2u + 3 > a hyperbola on the (y, H + 
plane. The case 4) p + V < 0, 2u; + 3 < is not realized as real solutions are absent. This 
result establishes that the intersection of the Friedmann 3-surface with the (fixed p, fixed x) 
2-plane is constituted in either one piece (ellipse) or two pieces (hyperbola). Thus in case 1) the 
allowed phase space is divided into two separate regions, the "upper" region where H + ^- > 
and the "lower" region where H + ^ < 0, and there is no way the trajectories can travel from 
one region to another. In case 2) these two regions meet along a 2-surface where H + = 0, 
and the trajectories can in principle cross from one region to another. In case 3) there are 
again two separate parts, now characterized by y > and y < 0, respectively. At first it may 
be difficult find a direct physical interpretation for the quantity H + ^- that characterizes the 
"upper" and "lower" region in cases 1) and 2), but it turns out that this combination is equal 
to the Hubble parameter in the Einstein frame, see Eq. (1391) . and thus the "upper" region 
corresponds to universes which expand in the Einstein frame, while the "lower" region has 
universes which contract in the Einstein frame. 

Related information can be also established by another approach. We may solve the Fried- 
mann constraint for H, Eq. ffTT]) . and then the condition for all variables to be real valued 
imposes an inequality 

(Mx)+ 3)^ + ^±ZW!>0. (2!) 
12x z ox 

In terms of physics this inequality can be interpreted as a restriction on the allowed values of 
y (see Table [TJ). There is no restriction in case 1), while the case 4) p + V < 0, 2u + 3 < 
is completely ruled out since no real solutions compatible with the Friedmann constraint exist 
in this case. Similarly, solving the Friedmann constraint for y leads to another inequality, 

°° 2 ((2uj(x) + 3)H 2 -^-^(p + V(x)))>0, (22) 



uj 2 (x) \ 3x 
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: fl^i, ff, p) = (2u(x) + 3) ^ + ^ w; = , (23) 



Table 1: Constraints from the Friedmann equation on the values of \I/ = y ( 12T1) and if (1221) . 

which can be interpreted as a restriction on the allowed values of H (given also in Tabled]). 
Analogously, once u(x) and V(x) are specified, we may get a third inequality from solving the 
Friedmann constraint for X clS well. 

In terms of the phase space geometry the inequality (1211) is saturated on a cylindrical 
3- surface 

y 2 k\ p + v (x)) 

12x 2 3x 

which is parallel to the H axis. In cases 2) and 3) it bounds the extent of the Friedmann 
surface (fl8l) in the x,y,p directions, i.e., when we project the Friedmann surface along the 
H direction to the (x,y,p) 3-plane, the projection lies within the bounds set by Bh- In case 
1) the projection covers the entire (x, y, p) 3-plane. Correspondingly, the inequality ff22l is 
saturated on the 3-surface 

B 4 : B^(x, y, H, p) ee {{2u{x) + 3)H 2 - ^^(p + V(x))^j = , (24) 

which is parallel to the y axis and bounds the extent of the Friedmann surface in the x, H, p 
directions. 

In cases 2) and 3) the Friedmann 3-surface T tangentially touches the cylindrical bounding 
3-surface Bh along a 2-surface defined by the union T U Br- The touching is tangential, 



as the trajectories can not go through Bh to the unphysical region. By substituting the 
condition fl23|) into the constraint equation (ITS!) we see that on the 2-surface T U Bh the 
quantity H + ^ vanishes, hence T U Bh joins the above defined "upper" and "lower" regions 
of the Friedmann surface. We may ask, again relevant in cases 2) and 3) only, whether the 



8 



trajectories do cross between these two regions, i.e. whether they do pass through JF U Bh- 
For this purpose let us introduce another 3-surface 

B H :B H (x,y,H,p) = H+^- = 0, (26) 

distinct from T and Bh, but with the property that its intersection with the Friedmann surface 
coincides with the union of T and Bh, i-e., T U Bh = T U Bh = Bh U From the scalar 
product 

,2 



K 



— ((l-w)p + 2V) (27) 



it follows that the trajectories pass through this intersection from "lower" region where H+^- < 
to the "upper" region where H + > if ((1 — w)p + 2V) |jRjg H > and in the reverse 
direction if ((1 — w)p + 2V) \^u6 H < 0- Only if w = 1 and V\ :Fu § H = do the trajectories, 
which are in the subspace T U Bh, remain there in their entirety, and thus block the passage 
between the "upper" and "lower" regions. 

In the next two sections we study the asymptotic behavior of solutions by finding the fixed 
points and their properties. Since the energy densities of different types of matter and the 
potential V(*&) scale differently under the change of the scale factor of the Universe, one is 
justified to consider different regimes separately, specified by the dominant component. 



4 Fixed points for potential domination (V ^ 0, p = 0) 

In the limit of vanishing matter density the phase space shrinks to three dimensions {x = 
ty,y = ^,H}, where the Friedmann constraint restrains the physical trajectories to lie on a 
two-dimensional surface T . We may solve the Friedmann constraint for H, as ( JTTT) . substitute 
it into Eq. (|7j), and thus in effect reduce the system 2-dimensional: 

( x = y 

{ v = {i - i) y 2 { + } + % 2 + y + iMr 3 (m*) - * f ) • 

This constitutes a projection of the trajectories on the original two-dimensional constraint 
surface in (x,y,H) to the (x,y) plane. The projection yields two "sheets": the "upper sheet" 
marked by the — sign, and the "lower sheet" marked by the the ,_|_, sign in Eq. ( 128|) . In 
three dimensions the former corresponds to the "upper" region where H + > 0, while the 
latter to the "lower" region where H + < 0. 

The consideration of only the variables {x, y} and Eqs. (128j) is in principle sufficient to 
follow the dynamics as the value of H can at each point be computed from the constraint 
f fTTj) . The only complication arises in cases 2) and 3) along the boundary y 2 = ^^ly+lf which 
marks the extent of the (x, y)-plane accessible for physical trajectories. In three dimensions 
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Table 2: Fixed points and their eigenvalues for the V ^ 0, p = case. 

this boundary corresponds to the border T U where + ^- = and the "upper" and 
"lower" regions meet. From the discussion around Eq. (I2T1) we know that with the exception 
of some special cases the trajectories do traverse from one region to another, and thus in the 
two-dimensional projection picture must change the "sheet" . However, along the boundary the 
Friedmann surface is positioned perpendicularly to the (x, y)-plane and the projection (12 8j) is 
then unable to encode the full dynamics of the system. Therefore with Eq. (128]) it is possible 
to follow a trajectory on one "sheet" forwards until it reaches the boundary or backwards until 
it starts from the boundary, while the step of changing the sheet remains to be accounted by 
relying on the full system (jSJ)-©- Still, the system ( 128]) is adequate for finding and describing 
the fixed points, at least as long as the prospective fixed points do not reside on the border of 
the two "sheets" . 

Standard procedure reveals that the dynamical system (l28l) is endowed with two fixed 
points, Table [2] lists their conditions and eigenvalues. The first fixed point satisfies 

^L.*.-2V(*.)=0, (29) 

which matches the second limit (|17l) . discussed in Sec. 2. It is rather surprising to note that 
a local extremum of the potential, jj^ = 0, does provide a fixed point only if the value of 
the potential vanishes at this point, V(^,) = 0. In the latter case the eigenvalues tell that 
if this externum is a maximum, < 0, the nature of this point is saddle, while a local 
minimum, ^p- > 0, is non-hyperbolic. But generally local minima with non-vanishing value of 
the potential do not figure as fixed points, while at the same time it is possible to have a fixed 
point residing on an arbitrary steep slope of the potential. This puzzling feature, however, 
becomes more meaningful when translated into the Einstein frame. As explained in Sec. 6 the 
condition ( |29l) gives a local extremum of the Einstein frame potential. 
The second fixed point satisfies 

1 1 dcu 







2w(**) + 3 ' (M#*) + 3 ) 2 d * 
i.e. exactly the same conditions (a)-(d) as the limit of general relativity for flat FLRW 
STT cosmology, discussed in the end of Sec. 2. Again, the eigenvalues are listed in Table 
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[2J In particular we see that on the "upper" sheet this point is an attractor if V > and 
(2tj+3) 2 ~ < 0' w hile on the "lower" sheet it can not be an attractor at all. (In 

principle it is also possible that the points and coincide. Yet, the properties of such a 
combined point are difficult to determine without knowing the exact form of u and V.) 

From Eq. ( TTTT) it is straightforward to compute that the values of H corresponding to the 
fixed points *. and ** are H, = ,j\ \J and H± = \J ^-^r~~, respectively. The 
result, which mimics de Sitter evolution in general relativity, was expected, since we saw in 
Sec. 2 that under the fixed point conditions ( J29l) and ( J30l) the full STT equations ©-(El) reduce 
to the equations of general relativity. 

So, having a model of STT with given u(^>) and V(^), one has to solve the conditions fl29|) 
and fl30|) to find the values of \l/ where fixed points occur. To determine the nature of these 
points it is necessary to compute the eigenvalues at these points, possibly using an appropriate 
limiting procedure. For the benefit of the reader, may we recall that when both eigenvalues are 
real and negative then we have a stable node (attractor), while real and positive eigenvalues 
indicate an unstable node (repeller). Complex eigenvalues with a negative real part give a 
stable focus (spiralling attractor), while a positive real part indicates an unstable focus. One 
positive and one negative real eigenvalue occur with a saddle point, but if the real part vanishes, 
the point is classified as non-hyperbolic and linear stability analysis is not enough to determine 
the behaviour of solutions around it. 

As an illustration let us consider an example 

3\? 

^ = 2(1-*) ' Vm = ~ 2{ * ~ °' 2)3 + 3( * ~ 0,2)2 ' (31) 
chosen to demonstrate some typical features that may crop up in a generic scalar-tensor theory. 
The shape of u and V, along with the phase portraits are shown on Fig. 1. To briefly go through 
the main features let us first recall that while the phase space is 3-dimensional (*, *, H) the 
physical trajectories lie on the 2-dimensional Friedmann surface which is comprised of two 
sheets. The domain \I/ G (0, 1] belongs to ) of Table [TJ meaning the values of \I/ are 

not limited, and the "upper" and "lower" sheets are separate. On the other hand, the domain 
* G [1, 1.7] belongs to the case 2), where the values of * are limited and the "upper" and 
"lower" sheet meet along the H + = line. Finally, the domain * e (1.7, oo) falls under 
the case 4) where the Friedmann constraint does not have any real solutions. 

There are three fixed points coming in two sets, one set for the "upper" and another for the 
"lower" sheet. Two of the fixed points satisfy the condition ( )29|) and happen to have the same 
corresponding nature on both sheets: the point at *, = 0.2 (where V = 0) is non-hyperbolic, 
while the other point at *. = 0.73 is a saddle. The third fixed point ** = 1 satisfies (130]) and 
is an attractor on the "upper" sheet, but a repeller on the "lower" sheet. Fig. 1 (up, right) 
illustrates how the trajectories in the * G [1, 1.7] domain start at the "lower" sheet repeller 
cross over to the "upper" sheet, and end up at the "upper" sheet attractor \Ev This behavior 



11 





Figure 1: Illustration to the example fl3T|) described in the text: the shape of to and V (up, 
left); 3- dimensional phase space with two sheets where the physical trajectories reside (up, 
right); 2-dimensional projection of the "lower" sheet (down, left), and "upper" sheet (down, 
right). 
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is reflected on the projection, which shows the "lower" sheet trajectories (Fig. 1, down, left) 
running to the boundary in order to just emerge on the "upper" sheet boundary (down, right). 



5 Fixed points for matter domination (V = 0, p ^ 0) 

In the case of cosmological matter (p > 0) and vanishing scalar potential the Friedmann con- 
straint restricts the solutions onto a three-dimensional surface in four phase space dimensions 
{x = ^$?,y = ty,H,p). However, the system is amenable to a change of the time variable, first 
used by Damour and Nordtvedt [12] in the Einstein frame, that allows to combine the field 
equations into a dynamical equation for the scalar field which does not manifestly contain 
the scale factor or matter density. In the Jordan frame this amounts to defining a new time 
variable p [15] 



dp = h c dt 



H + — 
2^ 



dt . (32) 



Then from Eqs. (J21)-(jll) the following "master" equation for the scalar field can be derived 
[ISIIM]: 

_ (M*) + 3)(l-w) 3 _ / 3(l + w) _ 1 ^2 

(+) 8^ 2 V 4 ^ 2u(^) + 3 dty J 

. 3(1 - w) . 3(1 -3w) , . 

+ — W - , \ r J —V = 0, 33 

(-) 2 (2cu($) + 3) ' V ; 

where primes denote the derivatives with respect to p. The transformation of the time coordi- 
nate (j32|) managed to align the phase space trajectories in a way that made a projection into 
two dimensions possible. Like in the previous section the upper signs in Eq. (133]) correspond 
to the "upper" sheet where the quantity H+-^ is positive, while the lower signs to the "lower" 
sheet where the quantity H + ^ is negative. 

The "master" equation can be relied on as long as h c = \H + ^\ is finite. At \I/ = the 
quantity h c diverges making the t-time to stop with respect to the p-time. Hence all t-time 
trajectories with finite ^ get mapped to Nl/' = 0, giving a false impression of a fixed point 
there. However, in Sec. 2 we concluded that ^ = comes with a space-time singularity and 
exclude it from present analysis. 

The other problematic points can be discussed by noticing that in terms of the new time 
variable p the Friedmann constraint (j5J) can be written as 

2 ft? p 

he — — r. (2^(io+3) v>2\ ■ (34) 



3$ 1 



To keep h c real, the right hand side of Eq. (|34|) must be nonnegative, thus constraining the 
dynamically allowed regions of the two-dimensional phase space (\1/, of the scalar field (see 
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Allowed range Allowed range 
of \i> of 



p > 2cj + 3 > < ^ 2 < oo 0<^ /2 <£g 



2) p > 2cj + 3 < 0<^ 2 < g| < < oo 



Table 3: Constraints from the Friedmann equation on the scalar field phase space in t-time 
from Eq. (12T1) and p-time from Eq. (134 



Table [3]). At the boundaries of the allowed regions the time transformation (|32|) fails to be 
meaningful. In case 1) the boundary 

12v|> 2 

2w(^) + 3 v ; 

is characterized by t-time stopping with respect to p-time, and as a result finite corresponds 
to \l/ 2 — > oo. Considering i-time to be physical permits us to exclude this boundary from 
analysis on physical grounds. In case 2) the limit ^' — > oo corresponds to finite since 
p-time stops with respect to i-time. As h c = \H + ^\ = the latter boundary also marks the 
joint of the "upper" and "lower" sheets. Therefore, in p-time it is not possible to follow the 
passage of trajectories from one sheet to another, but as the general considerations presented 
after Eq. (J27j) inform us, the trajectories do pass from the "lower" to "upper" sheet if w < 1 
and vice versa if w > 1, while in the w = 1 case the passage is blocked by trajectories which 
lie entirely on this boundary. Despite the shortcoming that the two-dimensional "master" 
equation ( 133]) is not able to capture all the details of the full dynamics in the four-dimensional 
phase space, we can still use it for finding the fixed points for \l/ as long as the fixed points do 
not reside on the problematic boundary. 

Let us consider the cosmological matter behaving like dust (w = 0) first. Introducing 
^ = x, ty' — z enables to write Eq. ( |33l) as a dynamical system: 

x' = z 

, + Ms)+3 „3 | 6u(x)+9-4x^ ^ 2 _ 3 3x (36) 

A ( — ) 8x 2 41 ~T 4a;(2w(x)+3) A ( + ) 2 4 2^(x)+3 " 

An argument completely analogous to the one put forth for Eq. ffT2|) . reveals a single fixed 
point, satisfying the conditions (a)-(d) dubbed as the limit of general relativity for flat FLRW 
STT in Sec. 2. The corresponding eigenvalues, to be evaluated at the fixed point coordinate, 
are given in Table HI In particular, this point is an attractor on the "upper" sheet if ^ > 0, 
while on the "lower" sheet attractor behavior is not possible. 

To complete the analysis it is important to verify that the fixed point just found in the 
p-time does indeed correspond to a fixed point in the cosmological t-time. But since we have 
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Case Fixed point Condition 



Eigenvalues 



w = * 
i 



1 

2o;(**)+3 



(2w(**)+3) 2 



7^0 



+ 

(-) 



32 



da; 



3 (2lj+3) 2 <W 



W 



none 



Table 4: Fixed points and their eigenvalues for the p 0, V = case. 

excluded the boundary (|35|) and consider only the trajectories with finite h c , it is immediate 
that = implies \? = 0, due to (|32|) . From Eq. (TlT!) now it also follows that at the fixed 
point the evolution of the universe obeys the usual Friedmann equation from general relativity, 
Hi = ^f-. This is expected, as the fixed point conditions were identical to the general relativity 
limit. 

The dynamical system in the radiation dominated regime (w = |) reads 
x' — z 

, + 2^)+3 3 | 8u^)+12-4x^> 2 _ (37) 

A ( — ) 12z 2 * 4x(2aj(a:)+3) ^ ( + ) 

There are no fixed points. For small values of \&' = z the system is ruled by friction on the 
"upper" sheet, as the — sign of the dominating term forces the vector flow to converge to the 
z' = axis. On the "lower" sheet, the the effect is the opposite (anti-friction). 



6 Discussion 



A lot of work in FLRW scalar-tensor cosmology has been performed in the Einstein frame 
which is obtained from the Jordan frame by two transformations [28| [30] : (1) a conformal 
transformation of the metric g^jy = ^g^fj, followed by a coordinate transformation to keep 
the FLRW form of the line element, dt = yftffdt, and (2) a redefinition of the scalar field 

(2w(tf) + 3)(d^) 2 = A^ 2 (d(f)) 2 . (38) 

The two frames are mathematically equivalent and hence also physically equivalent, as long as 
these two transformations are regular [281 129], El [30l EI] . Both transformations become singular 
in the limit ^ — > 0, while the latter transformation is also singular in the limit 2ui(ty) +3 — > 0. 
These singularities were scrutinized in Refs. [231 E2J E3] with the conclusion that both frames 
retain equivalence in the sense that in both frames the Cauchy problem fails to be well posed in 
these limits. In Sec. 2 we also observed in passing that approaching these two limits generically 
leads to a space-time singularity in the Jordan frame. 

The transformation (1381) happens to be singular in the limit of general relativity, 2^(^+3 ~" * 
0, as well. Ref. [20J also gives explicit examples of coupling functions where the phase 
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portraits for regions containing the GR limit are qualitatively inequivalent in the two frames. 
Therefore it is not guaranteed that the properties of the fixed point which resides at the 
general relativity limit are exactly the same in the Einstein frame. To actually compare the 
two frames in this respect, it would be necessary to perform a similar fixed point analysis in 
the Einstein frame. 

The other fixed point, however, should have the same properties in the Einstein frame 
as in the Jordan frame, as long as it resides in the region where the transformations between 
the two frames are regular. 

The conformal transformation relates the Hubble parameters in the two frames as 



He = ^= H j + — . (39) 




Therefore the sign of the quantity Hj + ^ which in our Jordan frame analysis distinguished the 
"upper" and "lower" sheet of the Friedmann surface, has a clear interpretation in the Einstein 
frame as indicating the expanding or collapsing universe. Also, as the scalar potentials in the 
two frames are related by 

VeM*)) = ^ W) , (40) 



it follows that the fixed point ^. given by Eq. f[29|) indeed corresponds to the local extremum 
of V E 



dV, 



E 



d<p 



1 



d^ 



d^ 

-j— =0. (41) 



9. 

Among the two types of fixed points we found, and it is the latter which deserves 
particular interest, since it satisfies the PPN limit of general relativity and good conformity 
with Solar System experiments is guaranteed. The existence of relies on the conditions 
(a)-(d) given in Sec. 2, while its nature (attractor or other) is determined by the eigenvalues 
listed in Tables |2] and HI Previous phase space analyses performed for specific examples of 
coupling functions and potentials which have considered this point [2H [20] are in accord with 
our general results. Still, we may also attempt a comparison with relevant studies carried out 
by other methods. 

Damour and Nordtvedt [12] argue that for large classes of coupling functions a 2 (0) = 
2w(l)+3 solutions are driven to the value of scalar field where the coupling function a vanishes, 
the scalar field decouples and we are left with general relativity. Their investigations were 
performed in the Einstein frame using approximate late time solutions and their results are in 
qualitative agreement with our results: for radiative matter, there is no specific fixed point, for 
dust matter general relativity is a fixed point. For instance taking linear coupling, a{4>) = Kef) 
with K = const > 0, in the matter dominated regime, they found that for K < | the system 
exhibits damped behavior, while for K > | the behavior is damped-oscillatory. In our analysis 
this can be compared to the ^ fixed point. Plugging in the coupling to the eigenvalues in 
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Table H] gives — | ± §(yl — indicating a stable node (attractor) for K < ^ and stable 

focus (spiralling attractor) for K > jq. The slight difference in the numerical factor may reflect 
the difference between the frames, or perhaps can be attributed to the fact that in making the 
approximation Damour and Nordtvedt drop some terms in the scalar field equation. 

Recently Barrow and Shaw [32] analysed the asymptotic behavior of homogeneous and 
isotropic cosmological solutions of general scalar-tensor theory containing a two-component 
perfect fluid: vacuum 'dark energy' (p = — p , essentially a constant potential) and sub- 
dominant matter density. They demonstrate that if STT cosmology evolves toward the de 
Sitter limit at late time, then there exists an asymptotic value of the scalar field \lZoo G (0, oo) 
such that for ^ — > it holds that to — > oo and — > for any e > 0. In our analysis of 
Sec. 4 this can be compared to the solutions closing in on the fixed point ^ yielding de Sitter 
expansion. The conditions for the existence of (a) - > and (c) ^+3)2 1^1*=** 7^ 
are slightly stronger than the ones of Barrow and Shaw. Namely, if ^ is identified with ^oo, 
then (a) and (c) imply the conditions given by Barrow and Shaw, but the reverse does not 
neccessarily follow, since (c) may not be satisfied. 

Finally let us note that scalar-tensor theory has some features which allow it to be seen 
as a simplistic toy model of the low energy effective actions derived from string/M-theory 
compactifications. Namely, the latter usually involve a number of scalar fields (moduli of 
the compactification) and some of these are coupled to the Ricci scalar in the Jordan/string 
frame description. For phenomenological reasons it is important to have the moduli stabilized 
at fixed values. Typical scenarios of string theory moduli stabilization involve generating a 
scalar potential whereby the moduli can stabilize at the minimum of the potential [35]. In 
our notation this parallels to invoking the fixed point However, it would be interesting 
to see whether the other fixed point, occuring at the singularity of the scalar field kinetic 
term, can also be generalized for string theory compactifications and what role can it play in 
stabilizing the string theory moduli. 

7 Conclusion 

We have considered flat FLRW cosmological models in general scalar-tensor theories with 
arbitrary coupling function ou(ty) and scalar potential V(^f) in the Jordan frame. Using the 
methods of dynamical systems we have described the general geometry of the phase space and 
found the scalar field fixed points in two distinct asymptotic regimes: potential domination 
(V ^ 0,p = 0), and matter domination (V = 0,p ^ 0). In nutshell there are two types of 
fixed points arising from different mechanisms: from a condition on the potential (equalling 
the local extremum of the Einstein frame potential) and from the singularity of the scalar 
field kinetic term. Approaching both types of fixed points the cosmological equations coincide 
with those of general relativity, yielding de Sitter expansion in the potential domination case 



17 



and Friedmann evolution in the matter domination case. However, for the Solar System 
experiments in the PPN framework only the fixed points of ^ type give predictions identical 
with those of general relativity. The nature of fixed points (attractor or otherwise) depends on 
the functional forms of u(^) and V(ty) according to corresponding eigenvalues given in Tables 
[2] and HI Therefore, in Jordan frame analysis, general relativity is an attractor for a large class 
of scalar-tensor models, but not for all. 

Provided the transformation relating the Jordan and Einstein frame is regular, there is 
an exact correspondence between the two frames and the Jordan frame phase space results 
should carry over to the Einstein frame. This is the case for the fixed point However, as 
the transformation of the scalar field fails to be regular in the limit of general relativity, the 
properties of the ^ fixed point may be altered in the Einstein frame. To establish whether 
or how the correspondence holds in this case calls for a separate matching investigation in the 
Einstein frame. 
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